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Abstract. We propose a method to implement a quantum memory for light based 
on ensembles of two-level atoms. Our protocol is based on controlled reversible 
inhomogeneous broadening (CRIB), where an external field first dephases the atomic 
polarization and thereby stores an incoming light pulse into collective states of the 
atomic ensemble, and later a reversal of the applied field leads to a rephasing of the 
atomic polarization and a reemission of the light. As opposed to previous proposals for 
CRIB based quantum memories we propose to only apply the broadening for a short 
period after most of the pulse has already been absorbed by the ensemble. We show 
that with this procedure there exist certain modes of the incoming light field which 
can be stored with an efficiency approaching 100% in the limit of high optical depth 
and long coherence time of the atoms. These results demonstrate that it is possible to 
operate an efficient quantum memory without any optical control fields. 
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1. Introduction 

Light is an ideal carrier of information both in the classical as well as in the quantum 
regime. To harness the full potential of light for quantum information processing it is, 
however, a major advantage to have access to quantum memories capable of storing the 
information from the light into an atomic medium and later releasing the information. 
Ideally such quantum memories should introduce as little disturbance as possible to the 
information encoded in the light field, i.e., the memory should be fully coherent and 
efficient. Furthermore it should be as simple to operate as possible. 

A large number of proposals have been developed for how one can construct 
quantum memories based on atomic ensembles [1]. As opposed to the most natural 
approach of storing one photon in one atom, the idea behind the ensembles based 
approach is to store photons into the collective states of an ensemble of atoms. Thereby 
one avoids the technical challenges associated with efficiently coupling single atoms and 
single photons. As a consequence the ensemble based approach considerably simplifies 
the experimental realization of quantum memories. A large class of quantum memory 
approaches based on atomic ensembles uses classical laser fields to control the memory 
process in such a way that the incoming field is mapped into the ground state coherence 
of the atoms. Based on this a number of key experimental advances have been achieved 
(See for instance [21 [3l HI [5l El [71 IHj ) . Here we shall pursue a different approach based on 
controlled reversible inhomogeneous broadening (CRIB) [9]. The original (transverse) 
CRIB is developed for impurities embedded in solid state systems. Due to the random 
orientations of the (static) dipole moments of the atoms, an external applied field will 
lead to different shifts in the transition frequency of different atoms. Hence an external 
applied field will dephase the polarization of the atoms. This will in essence turn off 
the reemission of the light absorbed into the ensemble which is therefore stored. The 
central idea in the CRIB approach is to later reverse the direction of the external field. 
If the atomic dipoles are fixed, the reversal of the field generates a shift in the opposite 
direction (see figure [He). After a certain time the polarization will rephase, causing 
an echo, where the light is reemitted. In the original approach an additional ground 
level was used and it was shown that the by driving this additional transition with laser 
fields one can obtain a readout in the backward direction which can in principle reach 
100% efficiency (within the theoretical models we will consider here, the efficiency is 
the only parameter characterizing the performance of the memory for a single incident 
mode [H [To], and we will therefore focus on this parameter below). From a practical 
perspective it is highly desirable to avoid the use of a third level and the associated 
lasers to drive that transition. 

Since the original CRIB proposal [9] (see also [11]) a number of modifications have 
appeared [121 [131 E] much experimental progress has been reported [151 lEl [iSl [13 
[T8l [T9l [20] . In particular a different longitudinal CRIB approach has been developed 
where a gradient of an external field is applied such that the shift depends on the position 
of the atoms [12]. For longitudinal CRIB it has been shown that one can construct an 
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efficient memory based on two level systems which in principle can reach an efficiency of 
100% [13] if the atoms have a sufficiently long coherence time. For transversal CRIB it 
was also shown that one can construct a quantum memory based on two-level systems, 
but in this case, however, the efficiency was found to be limited to 54% [21]. Here we 
develop the theory for an efficient quantum memory based on two- level atoms subject 
to transverse CRIB. We show that by varying the broadening in time so that it is not 
turned on initially but only applied during a short time interval as shown in figure [1] 
(d), the efficiency of the memory can in principle reach 100% if the optical depth of the 
ensemble is very high and the atoms have a very long coherence time. 

2. Setup and principle of operation 

Our system consists of two-level atoms that are coupled to a quantised electric field. 
We treat everything in one dimension and thus the electric field is only dependent on 
the spatial coordinate z. The atoms are assumed to be confined to a length L along the 
z-axis and are distributed with uniform density. The transverse extent of the ensemble 
and cross-section of the electric field mode are assumed to be the same with area A. See 
Figure [1] a) . 

The n'th atom has the two internal levels denoted \g)n and |e)„ for ground 
and excited state respectively. The two levels are connected by a transition with a 
matrix element p. We assume that the atomic transition frequencies are affected by 
both an intrinsic inhomogeneous broadening that we cannot control and a controlled 
inhomogeneous broadening that we can rapidly turn on and off and reverse in a time 
much shorter than the duration of the incoming pulse. These two types of broadening 
give rise to a detuning of the atomic transition from the incoming laser field of frequency 
cul. The two detunings for the n'th atom are denoted by Aon for the intrinsic broadening 
and A„ for the controlled broadening. Since we reverse the detuning the transition 
frequency is Ueg,n = ujl + ± A„ during the two periods where the broadening is 
applied, see figure [T]b). 

We are interested in storing light in two-level systems. Such storage is only possible 
if the lifetime of the excited state is very long. Therefore we shall assume that the excited 
state has a sufficiently long lifetime so that spontaneous emission can be neglected. The 
only source of decoherence of the atomic transition which we consider here comes from 
having an intrinsic inhomogeneous broadening of the atomic transition. Note that this 
assumption of negligible spontaneous emission from the excited state does not contradict 
the assumption that light couples to the transition. The quantum light we are interested 
in, couples to collective states which can have a much higher emission rates due to 
collective enhancement effects. The memory can thus have a large bandwidth even 
if the individual atoms have a negligible decay rate. Indeed as discussed below, the 
bandwidth of the memory is characterized by the quantity 
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Figure 1. Schematic representation of CRIB based memory operation, (a) Setup of 
the quantum memory. The incident Ught pulse is stored in the atomic ensemble of with 
length L and transverse area A and is later retrieved in the same direction, (b) The 
original atomic transition frequency of the n'th atom U!eg,n = '^l + Aon is only detuned 
from the electric field frequency wl by the intrinsic inhomogeneous broadening Aon- 
When the controlled broadening is turned on, the original excited state |e)„ is further 



shifted in energy either up (|e)n'^'') or down (|e)J^ '') so that LOeg^n = + Ao„ ± A„ 
(c,d) The width and direction of the controlled broadening is plotted as a function of 
time, (c) Standard approach where the broadening is applied long before the pulse 
enters the medium. A single reversal of the broadening initiates a rephasing which 
causes a photon echo to be emitted at a later time, (d) The procedure proposed here 
where the broadening is only applied in a short time interval of duration Td after the 
pulse has already entered the medium during the period Tp. After the storage time 
Tg the application of the opposite broadening for a duration Td initiates the retrieval 
process which last for a period Tp. 



where c is the speed of hght in vacuum and go = \/ huJi,/2ceoA. This bandwidth can be 
related to more experimentally accessible variables if we consider a specific model for 
the broadening. Throughout this work we shall assume a Gaussian distribution of the 
intrinsic broadening. In this case the decay of the collective atomic polarization due to 
dephasing is described by a Gaussian in time and we define the coherence time T2 as the 
1/e decay time of the atomic polarization. The precise definition and a more detailed 
discussion can be found in [Appendix A[ The memory bandwidth can then be related to 
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the optical depth do that will be observed when the controlled broadening is turned off, 
and the coherence time through = do/{T2y/7i). 

As it was shown in [221 ES], characterizes the fastest time scale for variations 
in the pulse shape that can be allowed if the memory is to store the light with high 
efficiency. For a memory based on three level A-type atoms controlled by a classical 
drive, any pulse shape with a pulse length Tp ^ l//i can be stored efficiently in the 
ensemble. Here we are, however, interested in two-level systems. It was also shown in 
[25] that there exist certain incoming pulse shapes of duration Tp ~ which can be 
efficiently absorbed and retrieved by having the light interacting only with two-level 
atoms. To control the storage and retrieval process it was proposed in [23] to rapidly 
apply a vr-pulse to transfer the population from the excited state to an auxiliary state 
for long term storage after the pulse was absorbed in the ensemble. Later another rapid 
71 pulse can be applied to initiate the retrieval process. It was shown that for the ideal 
pulse shape of the incoming field the efficiency of this memory protocol is identical to 
the highest obtainable efficiency for the memory. Here we shall follow a similar approach 
as the one suggested in [23] but without the use of a third auxiliary level. Instead we 
follow a CRIB-like approach [9], [21] and apply a broadening of the atomic levels. In our 
scheme the light pulse is first absorbed in the unbroadened ensemble. Afterwards an 
externally applied inhomogeneous broadening rapidly dephases the polarization of the 
atomic transition. This dephasing effectively turns off the reemission of the stored light 
allowing for its long term storage if there is negligible decoherence of the atomic state. 
At a later stage we assume that we can reverse the applied inhomogeneous broadening 
resulting in a rephasing of the atoms which effectively turns on the reemission of the 
pulse. As we shall show below, this process allows for an efficient quantum memory if 
the ensemble has a large optical depth. 

The memory protocol that we propose can be divided into 5 stages shown in 
figure Hid). In stage 1 the incoming pulse with duration Tp is incident on the atoms 
and gets absorbed without any controlled broadening present. In stage 2 with a (short) 
duration the controlled broadening is applied, which turns off the reemission of the 
excitation. We find that considerably higher efficiencies are obtained if we allow the 
incoming pulse to have a small tail into stage 2 and we thus allow this in the detailed 
simulations below. Stage 3 of duration Ts constitutes the main storage period. In this 
stage we assume that the external broadening is turned off. The storage duration Tg 
is only limited by the coherence time T2 of the atomic polarization. In principle this 
coherence time could potentially be prolonged by e.g. using spin echo protocols, but for 
simplicity we shall not consider this possibility here. Stage 4 marks the beginning of 
the read out process. Here the external broadening is reversed and applied for a time 
equal to the duration of the initial broadening in stage 2. The read out of the light pulse 
continues in stage 5 with duration Tp where the external broadening is again turned off. 

Most approaches to CRIB based memories [2T| [T2] do not include the stage 3 
introduced here but instead use longer stages 2 and 4. The idea behind this stage 
is to increase the flexibility of the readout. If the broadening is applied during the 
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whole storage period as in figure [T] a) the time it takes from we reverse the broadening 
to the pulse is retrieved is equal to the time from the storage until the reversal. Thus 
a long storage time gives a long response time of the memory. On the other hand, as 
soon as atoms are sufficiently dephased in stage 2, the reemission of the light is turned 
off and any further dephasing will not affect the performance of the memory. Hence we 
might as well turn off the broadening after the dephasing time which can be much 
smaller than the storage time Tg. Stage 3 cannot be longer than T2 if the memory is to 
be efficient but it can be arbitrarily short. If we want to retrieve the excitation stored 
in the atoms at any point in stage 3 we can turn on the reversed broadening (stage 4) 
and retrieve the excitations after a time Td which is independent of the storage time 
Tg. Additionally, small makes it easier to do numerical simulations of the memory as 
discussed below. 



In this section we present the model that we shall use to evaluate the performance of the 
memory. Here we shall be rather brief in describing the formalism. A more thorough 
derivation of the equations of motion and the principles used is given in P, [21]. To 
describe the light we consider the electric field operator given in the Schrodinger picture 



The atomic transition is subject to two types of broadening, the intrinsic broadening 
and the controlled reversible broadening. The probability distribution of the detunings 
is described by the functions Qq for the intrinsic broadening and Q for the controlled 
broadening. The distributions are assumed to be even functions and are normalized 



To derive the equations of motion we split the position and detunings of the intrinsic 
and controlled broadening into intervals of size zs, Aqs and respectively and define 
the atomic polarization operators 



Here the sum is over atoms having a position in the interval [kzs,{k + l)zs], an 
intrinsic detuning in the interval [/Aq^, (/ + l)Ao5], and a controlled detuning in 
the interval [mAo<5, (m + l)Ao<5]. The number of atoms that have positions and 
detunings in the mentioned intervals is N^im = {N/L)zsQq {l^os) ^osS (^^^A^) A^. The 
normalisation constant A/5 = {NLgop)/[hc) is here chosen such that all the constants 
of proportionality disappear in the final equations of motion below. In the limit 



3. Model 



by 

^(^)=^7o^siow(;2)e^"^^/^ + H.c. 



in terms of its slowly varying components that satisfy 



4iow(2;), 4\ow(^') =cd{z- z'). 



according to /_°;^ 6;o(Ao)dAo = 1 and /_°°^6;(A)dA = 1. 
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zs, Aqs, As — > 0, where <Skim,ij gets replaced by <Sij{z, Aq, A). The commutation relation 
is then 

4(z, Ao, A),S,A^', A'„ A')] = j^g^^^g^^Siz - z')S{A, - A',)S{A - A') 

X {Sji'Sij>{z, Ao, A) + 6j>iSi'j{z, Ao, A)) . 

We shall assume that all atoms are initially prepared in the ground state \g) and that 
only weak fields are incident such that most of the population remains in the ground 
state. When finding commutators we therefore approximate Sgg{z, Aq, A) ^ A5 and 
<See{z, Aq, A) ^ 0. Hcuce the commutator of Sge and Seg is a constant, and Sge can be 
regarded as proportional to an effective harmonic oscillator annihilation operator. 

If we introduce the slowly varying operator a{z, Aq, A) = e~^^^^^'^Sge{z, Aq, A) the 
Hamiltonian in the interaction picture and rotating wave approximation can be written 
as 

^^J^F r Iq (Ao + A) See{z, Aq, A) 



-9oP 



4iow(^)'T^(;2, Ao, A) + £l,^Jz)a{z, Aq, A)] } d^dAodA. 



This Hamiltonian describes the evolution during stage 2 where the controlled broadening 
is turned on. In stages 1, 3, and 5 where there is no controlled broadening we can find 
the Hamiltonian by removing the term with ASee{z, Aq, A) in the integrand. In stage 4, 
where the detuning is reversed we change the sign of this term. From the Hamiltonian 
we can find the equations of motion for the relevant operators. Assuming most of the 
population to be in the ground state the equations of motion are given by 

^ + c^^ 4iow(;^, t) = i^ 1^ P ^o(Ao)^(A)a(^, t, Aq, A)dAodA 
^^a{z, t, Aq, A) = -i(Ao + A)a{z, t, Aq, A) + i^^^^^8^,^^{z, t) 

Since these equations are linear and only couple effective annihilation operators we can 
ignore the hats on the operators and consider them to be equations of functions instead 
of operators. This allows us to calculate any normally ordered product of operators, 
and evaluate the efficiency of the memory. If we only consider a single incoming mode 
as we will do throughout this article, this efficiency is the only important parameter for 
characterizing the performance of the memory [1], [23] . 

Using the bandwidth of the memory /i defined above we can change to dimensionless 
units. We introduce the spatial coordinate z = z/ L and the time i = fi{t — z/c). We also 
introduce the dimensionless detunings Ao = Aq /fi and A = A/fi. Then we can define 
the dimensionless broadening distributions Go and G and the dimensionless electric field 
E. These are related to the old quantities by Go(Ao) = /i^o(/xAo), G{A) = /i^(/iA) and 
E{z,t) = £siow{z,t) / Using these dimensionless variables the equations of motion 
become 

Cj poo poo 

—E{z, i) = ij j G'o(Ao)G(A)a(z, i, Aq, A)dAodA, (1) 
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-j(t(z, t, Ao, A) = -i(Ao + A)a(5, t, Aq, A) + t). (2) 
at 

To describe reversal of the controlled broadening we can replace A by —A in ([1]) and 
dJ]). In the stages without any broadening (1, 3 and 5), we omit the term with A. 

To simplify the discussions below it will be convenient to introduce the polarization 
P for each value of the intrinsic detuning. This is defined by 

/oo 
G(A)a(5,t,Ao,A)dA. (3) 
-oo 

In stage 1 and 5 equations of motion ([1]) and (|2]) can be written in terms of P instead of 
a since the controlled broadening is not present (see (IB.ip and (]B.2|) in Appendix B). 



We also note that in the numerical simulations we shall assume that the broadening 
distributions are Gaussian, i.e. 

Go(Ao) = exp I ) , G(A) = ^1= exp ( ] . (4) 

Here the dimensionless widths are 7o = 7o//^ and 7 = 7/yU. 
4. Analytical theory 

Before looking at the results of the numerical simulations which provide a full assessment 
of the efficiency of the memory we can gain some intuition about the important aspects 
of the proposed protocol by doing simple perturbative calculations. Here we therefore 
investigate the process in stage 2 and 4 where the broadening first rapidly turns off and 
later turns on the absorption-reemission process in such a way that a large fraction of 
the excitation is left in the atoms during the storage interval (stage 3). To focus on this 
part of the dynamics we completely ignore the evolution in stages 1, 3 and 5 for now. 

For simplicity we shall neglect the intrinsic broadening so that our system is 
described by the equations 

^^E{~z,i) = i P G(A)a(5,t, A)dA, (5) 

^a(5,t, A) = -iKa{~z,t,K) + iE{~z,t). (6) 
ot 

We assume that all of the input pulse was mapped to P in stage 1 so that E{z = 0, t) = 0. 
At the beginning of stage 2 the initial condition is then (t(z, t = 0, A) = P(^)(5), where 
P^^^ denotes the polarization at the end of stage 1. To proceed we Laplace transform 
the equations [z — )■ u) and combine them into 

-^a{u, t, A) = -iA(T(M, t, A) - - / G{A')a{u, t, A')dA'. 
ot u 

We want to solve this equation perturbatively. To do this we introduce as{z,i,A) = 
e'^V(5, t, A) which satisfies the differential equation 

-.asiu, t, A) = — e'^* / G{A')as{u, i, A')e-'^''dA' . (7) 
at u j^^ 
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Under the assumption that as changes slowly as function of i we can replace asiu, i, A') 
in the integrand on the right hand side of ([7]) by as{u,i = 0, A') = P^^\u). Defining 
the Fourier transform of G by G(t) = G'(A')e~'^'MA' the approximate solution for 
a at the end of stage 2 is 

[u, A) = e-'^'^"'^ (^1 - ^ e'^*'G'(t')dt' j P^'^ (u) (8) 

For the case of Gaussian G we have G(t) = exp (— 7^P/2) describing the dephasing in 
time. 

The expression ([8]) gives the perturbative approximation for the evolution during 
stage 2 in the absence of an incoming pulse. In stage 4 we make the same approximation. 
The equations are the same as above but with A replaced by —A throughout to describe 
reversal of the broadening. Hence we actually look at the evolution of the polarization 
with reversed sign of A which evolves from the initial condition ([HD. At the end of stage 
4 the solution for a to first order in the perturbation is then 

C7^^\u, -A) = M - ^ y e'^*'G(t')dt' - ^ y e'^(*'-^^^)G(t' - /ird)dt'J P^'\u). 

From this expression and the definition we can find P. Using the assumption that 
G and hence G is an even function of its argument the total polarisation at the end of 
stage 4 can be written 

P^^\u) =(l-- (G{t')fdt'\ P^'\u). (9) 



u _ ^ 

The main idea behind the memory protocol we propose here is that excitations 
mapped into the memory during stage 1 are released during stage 5. The expression 
above describes how an initial excitation stored in the polarization of the atoms described 
by P^^'' at the end of stage 1 is mapped to the end of stage 4. To get an idea about the 
performance of the broadening mechanism we consider the efficiency with which this 
mapping is achieved. To do this we calculate the efficiency which can be expressed by 
1] = Jq \P^'^\z)\'^dz if we assume that P^^^ is normalized to unity (J \P^^\z)\'^dz = 1). 
Now we take the inverse Laplace transform {u z) of iQ. Assuming G to be Gaussian 
we can carry out the integration explicitly. Keeping only terms to lowest order in the 
perturbation we find 



r/(rd) = 1 - 2v^^erf(7rd) t P^^\z) f P^^\z')dz'dz. (10) 
7 Jo Jo 

This expression is the main result of this section. Ideally we would like to have an 

efficiency of unity. From the expression we see that an efficient memory operation can be 

achieved if 7 ^ /x, i.e., if the width of the broadening is much larger than the bandwidth 

of the memory. Alternatively this can be expressed in terms of the optical depth: as 

derived in [Appendix A| the measurable optical depth in the presence of broadening is 



d ~ v27r/i/7. Hence the limit where the broadening is efficient at storing the excitation, 
is equivalent to the limit where the optical depth is much smaller than unity after the 
broadening has been turned on. 
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The argument of the error function in f llUI) reflects the remaining reemission which 
has not been completely dephased at the end of stage 2 (the emission rate at the end 
is ~ drj/dTd). Hence if Td ^ I/7 the reemission is completely turned off and additional 
dephasing does not improve the efficiency considerably. Therefore there is no need to 
have the controlled broadening turned on during stage 3 as we said earlier. 

To have an idea about the validity of (11 01) . we can compare it to a numerical 



solution of ([5]) and ([6]) using the method described in Appendix B In figure [2] we plot 
the efficiencies calculated by the numerical and the perturbative approach. The two 
approaches are seen to be agree rather well for 7^/1. 

The results of [23] show that ideal performance of stages 1 and 5 is achievable 
for sufficiently dense samples. The results obtained in this section demonstrate that 
also stage 2 and 4 can be made to work, thus demonstrating that an efficient memory 
is achievable if we can apply a strong reversible broadening of the atomic transition 
7 ^ /i. To reach this result we have for simplicity assumed that the incoming field 
is only incident in stage 1. As mentioned above we, however, allow for an incoming 
field also in stage 2 of the protocol in our numerical evaluation of the efficiency of the 
protocol. As we shall see below, allowing for this small tail of the incoming pulse to leak 
into stage 2 improves the efficiency of the protocol and allows for a much more rapid 
convergence with the broadening 7 than predicted by the results of this section. 

5. Results 

Having verified that the memory can in principle work with near 100% efficiency in the 
ideal limit of dense samples and large broadening, we now turn to a detailed numerical 
investigation of the performance of the protocol for real parameters. Specifically we 
shall investigate the performance of the memory protocol for finite optical depth and 
finite broadening of the atomic transition. 

The details of the numerical procedure used to solve the equations of motion 



numerically are given in Appendix B We want to find the relation between the incoming 
and outgoing light fields. Since the underlying equations are all linear "beam-splitter 
equations" , which couple the annihilation operators of the effective harmonic oscillators, 
the relation between the incoming fields is also a beam splitter relation. As a result the 
connection between the incoming fields in stage 1 and 2 and the outgoing fields in stage 
4 and 5 can be written in the form 

Eoutii) = / Ke{1 t')EUf^Tn - t')dt', (11) 
Jo 

where tr = Tp + is the total duration of read-in sequence and also the duration of the 
read-out sequence. The detailed expression for the matrix kernel Ke can be found in 



Appendix B In the end we are interested in the efficiency of the whole process which is 



given by 



f-Mm 

V= I lEouM^di (12) 




Figure 2. EfSciencics of the broadening stages (stage 2 and stage 4) calculated 
numerically (dashed and dash-dotted lines) and perturbatively using (jlOp (solid line). 
In the figure we show the efficiency for an excitation stored at the end at stage 1 to 
still be present at the end of stage 4 as a function of the broadening 7 relative to the 
bandwidth of the memory ji. In all cases we assume that the stored polarization has 
the form P(^^(z) = 1 which represents the worst case scenario for the efficiency. In 
the figure we also compare two different durations of the broadening stages Td = l//x 
(dashed) and Td = 2//^ (dashed-dotted). The duration of the broadening stage has 
little influence on the efficiency as long as 7Td 3> 1. Hence in the perturbativc case we 
simply take the limit — > 00. 



if we assume that the incoming pulse is normahzed, i.e. /q^^^ \Ei^{t)\'^dt = 1. Using 
expression f lTT]) we can write the efficiency as 

r]= / E*^{fiTn-t')K,s{t',t")EUf^Tn-t")dt'dt", (13) 

Jo Jo 

where 

/•Mm 

J^eff(t',t") = / KUlt')KEii,t")di. (14) 
Jo 

In the numerical simulations we discretize the time and use a quadrature rule [23] to 
setup a matrix Keg{t,t') for the discrete times (nodes of the quadrature rule). With 
this matrix it is thus possible to obtain the memory efficiency for any vector containing 
the input fields at the discrete times. 

For evaluating the performance of the memory there are various approaches that 
one can take. From ( !T3|) we see that the memory efficiency depends on the shape 
of the incoming pulse. One approach to evaluating the performance is to look for 
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the incoming mode which has the highest efficiency. Since the final expression for 
the efficiency f|T3l) can be written in the form of a simple vector and matrix product, 
where the kernel matrix Kcs is self-adjoint, the maximal efficiency can be shown to 
be given by the largest eigenvalue of the matrix K^s. Alternatively the performance 
of the memory can be assessed by investigating how well the memory operates with a 
specially chosen mode that one may be interested in. Below we shall consider both 
these approaches. We emphasize, however, that there are many more methods of 
characterizing the performance. For instance the approaches that we take here do not 
characterize the ability of the memory to store multiple modes [25] and also ignore any 
information about the shape of the outgoing pulse, which may important for practical 
applications. 

To simulate the performance of the memory we consider a realistic situation where 
we have a collection of atoms with an intrinsic broadening of the optical transition 
described by the Gaussian distribution Gq defined in (jlj) with width 70 = 7o//i- Since 
we are interested in a quantum memory operating only on two-level atoms the width of 
this distribution inherently leads to a decay of the atomic polarization by an amount 
exp(— (rs/T2)^) with T2 = do/{y/T^fi) as shown in Appendix A At the same time, the 
width of the atomic line also determines the measurable optical depth c/q of the ensemble 
before the controlled broadening is turned on which is given by c/q = -\/27r/i/7o. In the 
investigations below we are mainly interested in how the performance scales with the 
optical depth of the atomic ensembles. We therefore fix the storage time to be 
equal to I/270 = T2/\/8 so that the maximal attainable storage efficiency is equal to 
exp(— rg^7Q) = exp(— 1/4) 0.78. With the storage time fixed relative to the dephasing 
time T2, the investigations of the dependence of the optical depth below essentially 
correspond to the scaling one obtains when varying the number of atoms in the ensemble 
while keeping all other parameters fixed. The allowed duration of the pulse to be stored 
is primarily determined by the duration Tp of stage 1. For a memory to make sense 
the duration of the pulse must be shorter than the memory time Tp < Tg. We therefore 
restrict ourselves to a duration Tp = Ts/4. Finally we chose a constant duration Td of 
stages 2 and 4. As shown above, the duration of this period is not so important as long as 
it is long enough that after stage 2, all the different polarizations are dephased sufficiently 
so that that negligible light gets out during stage 3. In the numerical simulations we 
are, however, constrained by having only a finite number of discrete frequencies. This 
means that there exists a finite time when the polarizations rephase again. Hence Td 
has to be chosen such that it is well below this rephasing time and we fix it at Td = l//i. 

In figure [3] we show the maximal efficiency obtainable for a given optical depth and 
width of the broadening. In the figure we see that the efficiency rapidly increases when 
we apply the broadening and saturates when the width of the applied broadening reaches 
a value around 7 > 3/i. Once the applied broadening reaches this value it is sufficiently 
broad to rapidly dephase the polarization and thus rapidly turn off the reemission of the 
absorbed light after the broadening has been applied. Furthermore we see that as we 
increase the optical depth of the ensemble do the efficiency approaches the maximally 
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allowed efficiency of 77 ~ 0.78. The procedure proposed here thus allows for efficient 
quantum memory operation using only two-level atoms. In particular by allowing the 
broadening to be turned on after the pulse is incident we are able to surpass the limit 
of ?7 = 0.54 identified in [21] for a two-level quantum memory based on the standard 
transverse CRIB approach where the broadening is turned on before the pulse enters 
the medium (figure [1] a) . Furthermore for a sample with a long coherence time T2 ^ r<j 
and a high optical depth ci ^ 1 we can in principle come arbitrarily close to an efficiency 
of 100%. 

0.8 I ^ ^ ^ ^ 1 




Figure 3. Efficiency rj for the optimal incoming mode and a Gaussian distribution of 
the controlled broadening. For a given set of parameters characterizing the memory 
we find the mode which gives the highest possible storage and retrieval efficiency and 
plot it as a function of the width of the controlled broadening "f/n for different optical 
depths do- Here do is the optical depth before the broadening is applied. Values of 
do in the legend are given in the same order as values of t] at j/fi ~ 10. In the 
simulation the storage time Tg is fixed at Ts = T2/V8, where T2 is the dephasing time. 
This limits the efficiency to 77 < exp(— 2(Ts/r2)^) w 0.78, which is shown as a dotted 
curve. For large optical depths do S> 1 and large applied broadening 7 > 3//i the 
curves approach the upper limit showing that we can have an efficient memory. The 
parameter ^ = do/{\/TTT2) characterizes the bandwidth of the memory and when the 
broadening is applied it reduces the optical depth to a value d ^ V27r/i/7. In the 
simulation we allow for an incoming pulse of duration Tp = Ts/4 and the dephasing is 
applied for a time ~ 1//^- 

Some examples of the optimal modes leading to the maximal efficiency in figure [3] 
are shown in figure HI For low values of the broadening, where the memory is inefficient, 
the optimal shape is a rather flat. As we begin to increase the width of the applied 
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broadening the optimal shape changes character. After a certain value of 7 which in 
this case is 7 ~ 1.4/i the shape of the modes for t < fiTp begins to resemble a Bessel 
like function, which is the shape of optimal modes identified in [23] (plotted in [26]) for 
the so called fast memory regime, which corresponds to our stage 1. This resemblance 
reflects that the mechanism in these approaches are highly similar. The only difference 
is that our dephasing mechanism which shuts off the reemission of the stored field is 
replaced in [23| by a vr-pulse taking the excitation from |e) to an auxiliary state. After 
the broadening is applied at t = /xTp the optimal mode shape rapidly drops to zero on 
a time scale set by the width of the broadening. This reflects that I/7 is the time scale 
needed dephase the excitation and thus terminating the read-in sequence. 
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Figure 4. The optimal mode shapes corresponding to the maximal efficiencies in 
figure [3] with do = 800. Here the controlled broadening is applied at t = /.jTp = 
do/(8\/27r) « 39.89, which is shown as a vertical dotted line, (a) optimal mode shape 
for different widths of the applied broadening. For 7 > 1.4^ the mode shape approaches 
the Bessel- like modes identified in [53] ■ (b) zoom in of (a) around /XTp. 

It is interesting to note that the efficiency in figure [3] approaches the stationary 
value much more rapidly when we increase the width 7//i than shown in figure [2J The 
reason for this difference is that we allow for a small tail of the pulse to leak into stage 2, 
which was not included in section HI In stage 2 the effective optical depth is given by 
d = y/27rfi/ -\/7^ + 7o ~ V^j^/l for 7 ^ 7o = V^j-J'/do. Hence for the regime where 
the memory is efficient 7 > /i the optical depth is below unity during stage 2. Regardless 
of this, the inclusion of a small optical field during this stage is still sufficient to alter 
the memory efficiency. This result emphasizes that the optical depth is not necessarily 
the correct physical parameter for characterizing the transient absorption of a pulse. 
The optical depth characterizes the fraction of the incident energy which is absorbed 
and not reemitted if the parameters of the memory are stationary. In our dynamic 
situation where the parameters are varying in time the optical depth does not correctly 
characterize the transient absorption of the pulse. In fact in our protocol most of the 
energy of the pulse is absorbed during stage 1 where the atomic line width is very narrow. 
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Hence the pulse that we store is much broader in frequency than the atomic hne width 
and most frequency components of the field see an optical depth much less than unity. 
This, however, only means that there is no stationary absorption and does not exclude 
that the field is absorbed and reemitted several times during the passage through the 
memory (e.g. the reduced group velocity of a pulse traveling through a transparent 
medium, such as glass, can be understood as a consequence of constant absorption and 
reemission events for the light. In this case the resulting pulse delay can be substantial 
even if the optical depth is very low). In this sense the idea behind the current memory 
protocol is to interrupt the frequent absorption and reemission events by applying the 
broadening which stops and later resumes this absorption and reemission process. 

Above we have focused on the maximal efficiency obtainable for the optimal mode 
shape. Since these optimal modes may not be available experimentally we can also 
consider how the memory performs for a given predefined mode shape. We therefore 
also calculate the efficiency for Gaussian input modes of the form 



This mode is characterized by two parameters, the center t^ and and the width and 
in figure |5] we show the efficiency as a function of these parameters. As we can see 
in the figure there is always a well defined maximum that we can find via numerical 
optimization. Furthermore in the limit where we expect the memory to work 7 ^ /i the 
best performance is achieved by sending in a rather narrow pulse of duration ~ A* 
right before broadening is applied. This results indicates that the best performance is 
achieved when the incident pulse resembles the sharply peaked optimal mode functions 
in figure HI 

If we optimize t^ and tc we can find the optimal Gaussian pulse for storage into the 
memory. The results of this optimization are shown in figure [6l Here we see that for the 
optimal Gaussian mode the efficiency also rapidly increases with the applied broadening 
and reaches a maximum for 7 ^ 3yU. Then efficiency actually starts to decrease which is 
a consequence of the Gaussian mode shape not resembling the changing optimal mode 
shape. The efficiency also increases with increasing optical depth of the ensemble do 
although much slower than for the optimal mode in figure [3l Unfortunately limited 
numerical resources prevent us from increasing the optical depth even further, but we 
believe that it will eventually converge towards the maximal attainable efficiency of 
V < exp(-2(rs/T2)2) ^ 0.78. Nevertheless our results are still able to surpass the limit 
of ?7 = 54% • 0.78 = 42% for the two-level memory protocol presented in [21] when we 
take into account the dephasing during the memory time. 

6. Conclusion and discussion 

We have proposed a method to make an efficient quantum memory for light based on two- 
level atoms. In our proposal the memory operation is controlled by applying and later 
reversing an external field which broadens the atomic transition. For sufficiently strong 




(15) 
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Figure 5. EfRciency of Gaussian modes for do ~ 800 and different widths tw and center 
tc of the Gaussian mode. The three different plots (a, b, and c) show the efficiencies 
for various different widths of the applied broadenings. Here tc = fiTp « 39.89 is the 
time at which the applied broadening is turned on (shown as vertical dotted line) . For 
7 3> the best performance is achieved for narrow pulses which are incident shortly 
before the broadening is applied. 

broadening and sufficiently high optical depth before the broadening is applied, we have 
shown that the memory operation can be limited only by the intrinsic decoherence (the 
T2 time) of the atomic transition. Most importantly, in contrast to most protocols 
the proposed memory can be efficient even without employing any optical control field. 
Since there is no need for an additional laser field we hope that this memory protocol 
may be simpler to implement in practice. 

Our protocol is not the ffist proposal for a quantum memory based on two-level 
atoms. In [T3| [T2] an efficient memory was proposed for the so-called longitudinal CRIB 
where a field gradient is applied such that the shift of the atomic transition depends 
on the position. On the contrary our protocol is based on transverse CRIB where the 
shift does not depend on the position. For this setting the only previous protocol [21] 
based on two- level systems had a maximal efficiency of 54%. For our protocol there is 
no such limit and the efficiency can approach 100% for sufficiently high optical depth. 
The efficiency is, however, dependent on the mode shape and converges to the ideal 
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Figure 6. Maximal efficiency of Gaussian modes (fTS)) for different optical depths 
before the broadening is applied do . Values of do in the legend are given in the same 
order as values of?7at7//z=10. The efficiency of a Gaussian pulse is optimized with 
respect to the width and center of the pulse and the figure shows the optimal efficiency. 
Similar to figure[3]the efficiency increases with the broadening and reaches the maximal 
value for 7 > 3fi. The efficiency also increases with increasing optical depth but the 
increase is much slower than for the optimal modes in figure |3l The dotted line at 
J] w 0.78 represents the maximal attainable efficiency due to the dephasing during the 
memory time. 



limit much more rapidly for the optimal mode shape than if we constrain it to be, e.g., 
a Gaussian. It would be interesting to study in more detail how the protocol developed 
here compares to the other memory protocols based on two level systems. The previous 
work, however, had a different scope than what we consider here and only studied the 
dependence on the optical depth after the broadening was applied and the effects of 
the intrinsic broadening were not included in the final results for the efficiency. Within 
our framework these previous studies thus correspond to the limit do ^ 00 and cannot 
directly be compared to our results. 

It has previously been argued that adding and reversing (transverse) broadening 
during the storage process generally reduces the memory efficiency for the optimal 
modes as compared to schemes which do not have this broadening [26]. In our scheme 
the broadening is absent during the period where most of the light is read into the 
memory. We therefore believe that for a given memory our proposed protocol will have 
the highest possible efficiency for the optimal mode. On the other hand it has also 
been shown that the previous CRIB protocols [211 [12] have a much larger multimode 
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capacity than protocol such as this one, where the broadening is not present during the 
storage [25]. This hmited muhimode capacity of our scheme is reflected in our results 
for Gaussian pulse shapes, where the efficiency becomes limited by the ability of the 
Gaussian to resemble the relatively few modes which are stored with high efficiency. 
In fact for the readout process we find that also for the case of an incoming Gaussian 
pulse, the outgoing mode shape resembles the outgoing mode shape for the optimal 
mode. Again this is most likely a consequence of the limited multimode efficiency: only 
a few storage modes are excited by the Gaussian, and these modes are later retrieved 
into outgoing modes which contain little information about the incoming pulse shape. 
This reshaping of the mode may be detrimental for some applications of the quantum 
memory, but for other applications the mode shape may be less important. For instance 
for quantum repeaters [27l |28] we are interested in interfering the output from two 
different quantum memories. Since both memories in this case will emit similar mode 
shapes, the reshaping of the mode is of minor consequence. For concrete applications 
of quantum memories a more detailed study will have to be performed to determine 
whether the increased efficiency of the present protocol outweighs the drawback of the 
limited multimode capacity. 
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Appendix A. Relation of the parameters of the model to physically 
measurable quantities 

In this section we relate the parameters of the model to physically measurable quantities. 
Specifically we shall express everything in terms of the optical depth in the absence of 
controlled broadening do and the coherence time of the polarization of the atoms T2. 
To do this we consider light propagating through the atomic ensemble without the 
controlled broadening being present. This situation can be described by fIB.ll) and 
( IB. 21) . If we formally integrate ( ]B.2|) under the initial condition P{z,t = — 00, Aq) = 
we get 




Taking the Fourier transform (t — > u) this equation becomes 




(A.l) 
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where PV reminds that we need to take the Cauchy principal value when integrating. 
Inserting ( lA.ip into the Fourier transform of fIB.ll) results in 

= ~E{z,uj) (^Go{uj) - FY Go(A[,)-^dA'o^ . 

From this we see that 

\E{2 = = exp(-27rGoM)|^(5 = 0,a;)|l 

From this expression we define the optical depth on resonance as do = 27rGo(0). If the 
controlled broadening is present this amounts to having a distribution of the sum of the 
controlled and intrinsic broadening Gsum instead of Gq. Hence we have d = 27rGsuin(0). 
If both Go and the controlled broadening distribution G are Gaussian with widths 7o//i 
and 7/yU respectively then Gsum is also Gaussian with width 7sum/yU = \/7o + 7^/^^- In 
the limit 7 ^ 70 we simply have 7sum ~ 7, so that 

dn = V2^-^ and d = v^^^^ ^ v^-^ 



70 V7o + 7' 7' 

To derive the temporal decay of the polarization we can assume that P in ([3]) is 
independent of Aq initially, i.e. P{z,i = 0, Aq) = Po(5). If no electric field is present, 
the evolution is given by (lB.2p and the total polarization becomes 

/oo poo 
Go(A)P(z,t, Ao)dA = NPoiz) / Go(A)e-^^«*dA. 
-00 J —00 

Hence for a Gaussian Gq given by @ we have 

-Ptotal(5, i) = Ptotal(5, t = 0) exp 

with T2 = V2/70 = do/ifiy/n). 



Appendix B. Details of the numerical method used to simulate the 
evolution. 

In this appendix we give the details of how we do the numerical simulations of the 
memory protocol. In short we Laplace transform the equations of motion in space 
{z u) and rewrite them such that we obtain an equation only involving the 
polarization. Then we discretize the broadening variables Aq and A such that the 
integrals become sums and convolutions become matrix products. The resulting vector 
equations have simple solutions expressed in terms of the matrix exponential. In the 
end we can find the Laplace transform of the electric field using the discretized and 
Laplace transformed versions of ([H) and ( ]B.1|) . Then the electric field in real space 
can be evaluated by applying numerical inverse Laplace transform which amounts to 
transforming the integration kernels. Detailed calculations for each stage are given 
below. 
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Appendix B.l. Stage 1 

In stage 1 the controlled broadening is not present so that a does not depend on A. 
Hence using the definition ([3]) we see that (T(2,t, Aq, A) = P(2, t, Aq). Also the term 
— iA(T(2, t, Aq, A) is not present on the right hand side of ([2]) so that we can write it 
together with ([T]) as 

= i P Go(Ao)P(5,t, Ao)dAo, (B.l) 

^P(5,t, Ao) = -iAoP(5,t, Ao) + i^(5,t). (B.2) 
ot 

Now we Laplace transform these equations in space (5 — )■ n). Combining the resulting 
equations and using the definition Ei^(t) = E{z = 0,i) we get 

d - ~ ~ 1/"°° i 

-jP(m, t, Ao) = -iAoP(M, t, Ao) - - / Go(A[,)P(m, t, A[,)dA[, + -EUt). (B.3) 

ot U U 

We choose discrete values of Aq that we call Aqi, . . . , Aoi^-. For simplicity we 
choose them such that they have a constant step A05. After discretization (IB.SP becomes 

d ~ ~ ~ I ^ ~ i 

-^P(u, t, Ao,) = -iAo,P(n, t, Aq,) - - V Go(Aofc) Ao5P(u, t, Aqa.) + -^in(t) (B.4) 

fe=i 

To simplify the notation we wish to write the discretized equations in vector form. 
To do this we define vectors Aq, V{u,t) and go with elements given by (Aq), = Aoj, 
(P(M,t))j = P(M,t, Aoj) and (go)j = Aq^GoI^oj) for 1 < j < -f^- Also for any vector v 
we define a matrix -D(v) with elements D{\')jj = vj for all j and D{\)jk = for j 7^ k, 
i.e. a diagonal matrix with v as its diagonal. If we further let h*^^^ be a i^-dimensional 
vector with constant elements hj = 1 for all j then we can write flB.4p as a vector 
equation 

^P{u,i) = -iD(Ao)P(n,t) - -h(^^)go^P(n,t) + -i?in(t)h(^^ (B.5) 
ot u u 

Note that h'^-^^gj^ is a matrix with each row equal to the vector gg. Defining 

Mi(n) = -iD(Ao) - -hWgo^ 

n 

we can write (IB.SP as 

4p(w,t) = Mi(«)P(M,t) + -Ein(t)h(^). (B.6) 

In stage 1 the initial condition is that the memory is empty, i.e P(5, t = 0) = 
or equivalently P(M,t = 0) = 0. The solution to flB.6P can then be expressed using the 
matrix exponential 

ft 1 

P{u,i)=i -eMMi(u)t')h^^^Ei,,{i-t')dt'. (B.7) 
Jo u 
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Appendix B.2. Stage 2 

The memory protocol works by having most of the light absorbed during stage 1. 
Nevertheless we still allow for some light also during stage 2. Hence stage 2 is described 
by ([I]) and ([2]) with the initial conditions 

E{z = 0,i) = E;^{i+ fiTp), 

a(z,t = 0,Ao,A) = a«(5,Ao,A) = pW(z,Ao) 

where a^^^ and P^^^ are two ways to denote the polarisation at the end of stage 1; and 
Ei^ is the same function that was defined for stage 1. Here and in the beginning of all 
the subsequent stages the time t represents the time since the beginning of the current 
stage. This resetting of time is convenient since the equations of motion change from 
one stage to another. Since we want to consider Ei^ as being an input pulse for both 
stage 1 and stage 2 we have to shift its argument by the duration of stage 1 when we 
solve the equations of motion in stage 2. 

Taking the Laplace transform of ([I]) and (I2l) and combining the resulting equations 
as for stage 1 we get 
d ^ ~ ~ 

-j(T(M,t,Ao,A) = -i(Ao + A)a(M,t,Ao,A) 
ot 



1 poo poo 

- / Go(A[,)G(A')aKt,A[„A)dA[,dA' (B.8) 

^ J —OO J —OO 

i 



U 

Now we have to split both detunings Aq and A into respectively K discrete values 
Aoi , . . . , Ao_ft: and discrete values Ai,...,Aiv. Again we assume that they have 
constant steps Aqs and A^ respectively. Defining aj^k{u,i) = a{u,i, Aqj, A^), the 
discretized version of fIB.SI) can be written as 

d ~ ~ ~ 

-^aj±{u,t) = -i(Aoj + Afc)crj-fc(n,t) 

^ K N 

~-J2Yl Go{Aoj,)G{Ak')AosAsaj,,kiu, i) (B.9) 
j'=l k'=l 

u 

We see that we can write this expression in vector notation as before if we combine 
the indices j and k into one. Hence we define vectors Aj., (T{u,t) and g with 
elements (A±)(j_i)7v+fc = Aqj ± A^., {a{u,i))^j_i)M+k = ^j,k{^^^) and {g)^j_i)N+k = 
AosAsGo{Aoj)G{Ak) for 1 < j < A' , 1 < A; < A^. With these definitions (iROj) becomes 

^CT{u,i) = -iDiA+)a{u,i) - ih(^^)gV(n,t) + + /irp)h(^^). (B.IO) 

at u u 

Defining M2{u) = -iD{A+) - h^^^^g^/w we can write (iRTOj) as 

= M2iu)cT{u,i) + -Ei„(t + /irp)h(^'^) (B.ll) 
ot u 
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and the solution at the end of stage 2 is given by 

= - / exp(M2(n)t)h(^^)EiJ/irR - t)dt + exp(M2(u)fiTd)(T^^\u) (B.12) 
u Jo 

with tr = Tp + Td as defined earher. Here we have used the initial condition cr^^^u) which 
in the vector notation has elements {cr^^\u))(j-i)N+k = (Piu, fiTp))j ior 1 < j < K and 
1 < A; < with P(m, i) given by (fRT]) . 

Appendix B.3. Stage 3 

To describe stage 3 we could in principle have used the same equations as for stage 1 but 
then we would have lost the information about how the individual frequency components 
of the controlled broadening evolved. Hence doing this would not have allowed us to 
describe the rephasing in stage 4. Instead we use an equation similar to (IB.9P and set 
Afc = 0. Here we do not have an incoming field so that the evolution is given by 

—^ajk{u,t) = -iAojajk{u,t) ^^Go{Aoj')G{A'^)Aos^so-j'k'{u,t). (B.13) 

^ j'=i k'=i 

We define a vector A with elements (A)(j_i)7v+fc = Aq-,- for 1 < j < i^, 1 < A; < and 
the corresponding matrix M^(u) = —iD(A.) — h.^^^^^ /u. Then we can write ( IB.lSp as 



^(t{u, t) = M^{u)ct{u, t) (B.14) 
ot 

so that the solution at the end of stage 3 is given by 

(T^^\u) = exp{Ms{u)iJTd)(T^^\u). (B.15) 

Appendix B.4- Stage 4 

For stage 4 we define M^i^u) = — iZ}(A_) — h^^^^'g^/u and the evolution of a is then 
described by an equation of the same form as flB.14p with the solution for a given by 

cr{u,i) =exp(M4(n)t)cr(3)(n). (B.16) 

From this expression we can find the electric field using the Laplace transformed and 
discretized version of ([1]) which can be written as 

E{u,i) = -g^(T{u,i). (B.17) 

u 

We define the matrix 

jiKNxKN)^^^ = exp(M3(M)/irs)exp(M2(M)/ird) 
and the kernel 

h{u,i,t') = ^g^exp(M4(M)t)exp(M3(M)/irs)exp(M2(M)t')h^^^^ 

Then we can use f lRT2|) . flRTSj) . f lRT6|) and f lRTT]) to write the field as 



E{u, t) = / hiu, I t')EUf^T^ - t')dt' 

JO (B.18) 
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Here the first term describes how the field incident in stage 2 is read out during stage 4 
and the last term describes the readout of the field incident during stage 1. 

Before we can write the second term in this expression as something involving the 
input field we need to express cr^^) (u) in terms of P*^^-* (u) so that we can use (IB. 70 . In the 
vector notation P^^\u) is a ii'-dimensional vector while cr^^^u) is a i^A^-dimensional 



vector as defined in the end of Appendix B.2 This definition means that every element 



in block number j of cr^^^u) with length has the same value (P^^\u))j. Hence when 
the matrix J^^^^^^\u) is multiplied with (t^^\u) each element of the resulting vector 
is given by 

KN 
k=l 

K N 

j'=l k'=l 

In the last line (cr(^)(n))(j/_i)Ar+fc/ = {P^^\u))ji so that it can be taken out of the k' 
summation. Hence if we define a matrix J^^^^^\u) with elements 

N 

k'=i 

then J'^^^''^^\u)(t'^^\u) = j'^^^''^\u)P^^\u). Then if we introduce another kernel 

Uu,i,t') = -^g^exp(MJu)i)J^'^''''''\u)exp(MAu)t')h^''^ 

and take the inverse Laplace transform (u — )■ 5 = 1) we can finally express the output 
field during stage 4 in terms of the input field 

Eont{t) = / h{i, t')Ei^{fiTn - t')dt' + / k2{i, t')E,^{fiTp - t')dt' (B.19) 
Jo Jo 

where the kernels ki and ^2 are inverse Laplace transforms of ki and ^2 evaluated at 

z = l. 

Appendix B.5. Stage 5 

In stage 5 the initial condition for a is given by f lB.16P with t = fiT^- Defining 
j^^iKNxKN)^^^ = exp(M4(w)^rd)exp(M3(n)/ir,), 
^(KNxKN)^^^ = exp(M4(M)/ird)exp(M3(M)/irs)exp(M2(M)/ird), 



we can write the initial condition as 

j^iKNxKN)^^^ exp(M2(M)t')h('^'^)^in(/irR - t')dt' 



u 



The electric field can be found using the Laplace transformed and discretized version of 



E{u,t) = ^g;^p(u,t), 
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where P{u,t) = exp{Mi{u)t)P^'^\u) is the solution to flR6i) with = 0. The initial 
condition P^^^(u) is defined by ([3]). In vector notation this translates to blockwise 
summation of a^'^\u) with respect to the index of the controlled broadening weighted 
by the appropriate value of the controlled broadening distribution. We can apply this 
blockwise summation to the matrices L^^^^^^\u) and B^^^^^'^\u) directly. Hence 
we define L^^^'^^^u) and B^^''^\u) with 

N 
N N 

j'=i k'=i 

where we also apply the same reduction for the matrix B^^^^^^\u) as for j(^^xX7V)^^) 
previously. We can now use kernels 

h{uXt') = --^g^exp{M,{u)i)L^'''<''''\u)exp{M,{u)t')h^''''\ 
hiu,i,t') = -^g^exp(Mi(«)t)5(^x/^)/ ) exp(Mi(M)t')h('^^. 

to express the output electric field in terms of the input field Ei^. Taking the inverse 
Laplace transform (m — >■ z = 1) the output field is 

Eontii) = / h{i, t')EUfiTn - t')dt' + / t')E,^{fiT^ - t')dt'. (B.20) 

Jo Jo 

where the kernels and are inverse Laplace transforms of k^ and k^ evaluated at 

z = l. 



Appendix B.6. Concluding remarks 



With the above results we have found the sought relation between the input and output 
field. The final expressions in fIB.lQP and (IB.20p give the output field during stages 4 
and 5 respectively. They can be combined into a single expression f|TT]) if we define 



KEii.t') 



i < flTd,t' < flTd] 
i < /iTd, t' > fiTd] 

i > fiTd,t' < fiT^; 

i > fiT^, t' > fiTd- 



ki{t,t') for 

k2(t,t' - fiTd) for 

k^{i-^TA,t') for 

k^it — /iTd, t' — /iTd) for 

Here the shifts of the arguments of the kernels are performed because we go from 
considering E^nt to be output fields only in stage 4 or only in stage 5 (where time was 
reset at the beginning of each stage) to having a single function E^nt giving outgoing 
field in both stages. In actual simulations we use expressions (1B.19P and (lB.20p directly, 
but expression f|TT]) is a more convenient formulation for the arguments made in the 
main text. 

For inversion of the Laplace transform we use numerical integration along the Talbot 
contour [29]. Since this method has exponential convergence, only a few discrete Laplace 
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moments u are needed. This makes it possible to compute and store certain intermediate 

results for each u. For example the kernels ki, ^2, and ki can be computed much 

faster if the matrices Mi(-u), M2{u) and M4^{u) are diagonalized. Hence we compute and 

store the eigenvalues and eigenvectors of those matrices, so that matrix exponentials 

can be evaluated by exponentiating the eigenvalues instead of the whole matrix. 
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